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The accuracy of the procedure that converts the experimentally determined electron spectrum 
, associated with fission of the nuclear fuels 235 U, 239 Pu, 241 Pu, and 238 U into the v e spectrum is 

• examined. By using calculated sets of mutually consistent spectra it is shown that the conversion 

procedure can result in a small ~1% error provided several conditions are met. Chief among them 
t are the requirements that the average nuclear charge (Z) as a function of the (3 decay endpoint 

. _ ■ energy is independently known and that the v e spectrum is binned into bins that are several times 

jij ' larger than the width of the slices used to fit the electron spectrum. 

< 

^ ; I. INTRODUCTION 

Good knowledge of the v e spectra of nuclear power reactors is an important ingredient of the study of 
neutrino oscillations with reactors The next generation of oscillation experiments (e.g. Daya-Bay 
and Double Chooz aim at an unprecendented accuracy in their search for the mixing angle #13. Even 
though the multi-detector scheme used in them will greatly reduce the dependence on knowledge of the 
reactor spectrum, a better understanding of it, and the corresponding error bars, would be clearly beneficial. 
Nuclear reactor monitoring with neutrinos is another potential application of the reactor v e detection 0]. 
Again, good understanding of the spectrum and its errors is a necessary condition for its success. 

The power in nuclear reactors originates in the (time dependent) contribution of fission of four nuclear 
fuels: 235 U, 239 Pu, 241 Pu, and 238 U. During a typical fuel cycle the dominant contribution of 235 U at the 
\ beginning of the cycle slowly decreases, and the contribution of the reactor produced nuclei 239 Pu and 241 Pu 

. increases. While 238 U represents ~ 97% of the fresh fuel rods, its fission, caused only by the fast neutrons, 

contributes only about 10% of the reactor power and changes little during the refueling cycle. In order to 
qq , determine the reactor v e spectrum, and its time development, one has to know, therefore, the corresponding 

v e spectra associated with the delayed (3 decay of the neutron rich fission fragments of these four fuels. 
f**» , There are two principal approaches to the determination of the v e spectra of the individual nuclear fuels. 

One of them combines knowledge of the fission yields, i.e., the probabilities of populating the individual 
fission fragments characterized by their mass A and charge Z, with the knowledge of their, often complex, 
(3 decay. The D e spectrum is then the sum of the contributions of all (3 branches of all fission fragments. 
The main source of uncertainty of this method is the fact that the (3 decay characteristics of some fission 
fragments are unknown (or poorly determined), particularly for those with very short lifetimes and hence 
high Q values. It is then necessary to use various nuclear models in order to characterize the (3 decay of 
these "unknown" nuclei. 

The other approach is based on the experimentally determined electron spectrum associated with fission 
of the individual fuels which is "converted" into the v e spectrum. The resulting spectrum uncertainty is 
then a combination of the experimental uncertainties of the electron spectrum and the uncertainty of the 
conversion procedure. The accuracy of the conversion procedure is the topic of the present work. 

The electron spectra corresponding to the thermal neutron fission of 235 U, 239 Pu and 2 Pu were measured 
in a series of experiments using the beta spectrometer BILL at the ILL High Flux Reactor in Grenoble [f| 
in 1982-1989. In the quoted papers the electron spectra were converted into the v e spectra. These spectra, 
supplanted by the calculated spectrum of 238 U (typically from Ref. Q ) , were used in the analysis of essentially 
all oscillation experiments performed so far. At short distances from the reactor they agree quite well with 
the measured v e signal [l|. 
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The conversion procedure used in Refs. [5[ is only briefly described there. The corresponding error is 
estimated to be 3-4% with a weak energy dependence. Since the statistical and normalization errors of the 
electron spectra usually exceeded the estimated conversion uncertainty, its precise value was less important 
so far. Here, I wish to return to this issue, two decades later, in order to determine the uncertainty in the 
conversion procedure more firmly, and study its possible optimization. 



II. CONVERSION FOR ONE NUCLEUS 



Throughout this work it is assumed that all relevant /3 decays have the allowed shape. The possible effects 
associated with the unique first forbidden decays (and higher order forbidden decays) are neglected. Thus, 
the electron spectrum is assumed to be of the shape 

N(E e ) = k(E , Z)E ePe (E - E e fF{Z + 1, E e ) , (1) 

where k(Eo, Z) is the normalization constant, E e and p e are the full electron energy and momentum, and 
F(Z + 1, E e ) describes the Coulomb effect on the outgoing electron. The v e spectrum as a function of the 
antincutrino energy E v is obtained from eq.([T]) by substituting E v = Eq — E e in the above formula. 




FIG. 1: (Color online) Example of complex electron and corresponding i/ e spectra. This example is for a hypothetical 
/3 decay of a Z — 45 nucleus with 40 branches. The endpoints are randomly distributed, with average spacing of 200 
keV, and the largest kinetic energy of the electrons is 8 MeV. The branching ratios are also randomized. 



For a single branch decay, therefore, the conversion procedure is a trivial one. Note the basic differences 
between the electron and P e spectrum in that case. There is a finite probability of emission of an electron 
with vanishing momentum. Thus, the electron spectrum begins at a finite value at E e = m e (p e — 0), while 
the v e spectrum vanishes for E v = 0. On the other hand, near the endpoint the electron spectrum vanishes 
as E e — > Eq, while the 9 e spectrum approaches a finite value as E v —> Eq — m e . Moreover, since electrons 
are attracted to the positively charged nucleus the electron spectrum reaches its maximum at lower energies 
than the v e spectrum. 
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An illustrative example of a complex spectrum is shown in Figure [T] Note that the full electron energy 
is used on the x-axis; hence the electron spectrum begins at m e . Due to the differences listed above the 
smooth electron spectrum is associated with a jagged v e spectrum. That is most noticeable at low energies, 
while at higher energies the two spectra are quite similar. (The very upper edge is not shown. Naturally the 
electron spectrum extends m e further than the v e spectrum.) 




FIG. 2: (Color online) The fit to the complex spectra in Figure Q] In the upper panel the smooth electron spectrum 
is fitted. The size of the slices and the corresponding notation is indicated. The plotted quantity is the normalized 
deviation of the fitted spectrum from the original one. In the lower panel the endpoints and intensities of the fitted 
branches from the electron spectrum are used to reconstruct the spectrum. Its normalized deviation from the true 
original i/ e spectrum is plotted for two slice intervals (300 keV - full line and 500 keV - dashed line). Note the change 
of scale between the panels. 



Let assume that the electron spectrum in Fig. [T]has been experimentally determined, but the endpoints 
and intensities of the individual branches remain unknown. How can one convert the electron spectrum into 
the v e one? The procedure used in Refs. [f| is sensible and practical. The electron spectrum is divided 
into a number n of slices. In Q these slices are all of equal width, and that choice again is simplest and 
economical, since it minimizes the number of free parameters. Then, starting with the highest energy slice, 
one assumes the allowed shape, Eq.([T]), and fits for the corresponding cndpoint and intensity. Once these 
parameters were determined, the resulting branch spectrum is extended down to E e = m e and subtracted 
from the full spectrum. The procedure is then repeated for the next slice, etc., until all n slices (or assumed 
branches) are characterized by their endpoints and intensities. The free parameter n, the number of assumed 
branches, is then optimized so that the true original experimental electron spectrum, and the fitted one, are 
as close to each other as possible. 

Since the electron spectrum is always a smooth function of energy, the procedure works well, as demon- 
strated in the upper panel of Fig. [5] For slices of 0.3 or 0.5 MeV width the maximum deviation is less than 
1%, except for the lowest and highest electron energies. (Since the spectrum is assumed to be measured in a 
finite number of points, every 50 keV in this example, the number of slices n is constrained also from above.) 
Having approximated the complex f3 decay scheme by n fictitious branches characterized by their endpoints 
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and branching ratios, one can construct the corresponding v e spectrum. Comparing it to the true original 
spectrum, one can see that the fit is less perfect (see the lower panel in Fig. [2]) . That is not surprising. The 
steps associated with true endpoints cannot be faithfully reproduced by the fit, hence the jagged form in 
Fig. O For practical application this is less important; if the v e spectrum is binned in bins wider than the 
slices of the fit, the up and down features in Fig. [2] would be smoothed out. 

III. SPECTRUM OF FISSION FRAGMENTS 

The fission fragments are broadly distributed in nuclear masses and charges. For example, the calculations 
used in Ref. @ were based on 1153 different fission fragments ranging from Z, A = 23,66 to 72,172. Most of 
them are radioactive, and many have complex (3 decay schemes. The resulting electron and v e spectra are 
therefore superpositions of many thousands of (3 decay transitions. 
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Electron full energy and antineutrino energy (MeV) 



FIG. 3: (Color online) Calculated spectra of electrons and u e corresponding to the thermal neutron fission of 235 U at 
equilibrium. 

An example of the calculated electron and D e spectra is shown in Fig. [3] They correspond to the thermal 
neutron fission of 235 U using the data files in Ref. @. Clearly, within the rather coarse scale of the figure, 
the two spectra are essentially identical; thus the conversion procedure should work well. Both spectra, 
calculated with a fine step of 10 keV, appear to be smooth. In the case of v e this is so because so many /? 
decays contribute that the step at the endpoint of each of them is not noticeable. 

The equilibrium spectrum in Fig. [3] is calculated using 

= Z A .zY(Z, A)XAE ePe {E l - E e ) 2 F(Z + 1, E e ) , (2) 

where Y(Z, A) are the cummulative fission yields, and arc the (3 decay branching ratios of the ith branch of 
the fragment Z, A with endpoints Eq . The v e spectrum is obtained from the same formula by the substitution 
E v = Eq — E e for every branch. 
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As stated above, the allowed shape is assumed for all j3 branches. This assumption should not affect the 
conclusions of this work, which is based on the consistent comparison of the converted v e spectrum with 
the calculated spectrum using Eq. (J2J). On the other hand, the actual spectrum associated with fission will 
clearly contain a number of forbidden (3 transitions. Among them the nonunique first forbidden ones have 
usually allowed shapes anyway, and the higher order forbidden decays could be safely neglected. It is difficult 
to estimate the error associated with the neglect of the known shape factor of the unique first forbidded 
decays, but it is not expected to be significant, since the converted v e spectra agree well with measurements 
at close distances from the reactor 

Now, let us use the conversion procedure described above and convert the electron spectrum, assumed 
to be known accurately, into the v e spectrum. When the smooth electron spectrum is divided into n slices 
of sufficiently narrow width, each slice can be well reproduced as a piece of the allowed spectrum, Eq.fl]). 
There is only one parameter to fit for each slice, the endpoint Eq, while the normalization is determined 
automatically. The dependence on the nuclear charge Z is rather weak, a good fit is obtained with a constant 
Z. 




Antineutrino and electron full energy (MeV) 



FIG. 4: (Color online) Deviations of the fit for the electron (dashed line; deviations very small and not visible in the 
plot) and Ve spectra, calculated with a constant Z = 47. The fit with 100 keV slices is the jagged line; interpolating 
the antineutrino spectrum over 250 keV wide bins the smoother line results. The spectra are for the thermal neutron 
fission of 235 TJ at equilibrium as in Fig. [3] 

Having the set of n endpoints and branching ratios, one can again reconstruct the v e spectrum and compare 
it with the v e spectrum derived from eq. |T]). There one encounters two problems. First, instead of many 
thousands of branches, there is now a much smaller number n of branches. Hence, the steps at the endpoint 
of each branch become more pronounced in comparison with the true v e spectrum. Second, the choice of the 
nuclear charge Z becomes important, unlike in the fit of the electron spectrum. These features are illustrated 
in FigHJ There, the fit was performed with a constant Z = 92/2+1 = 47, with slices 100 keV wide which 
were also rebinned to 250 keV wide bins. The electron spectrum, with or without binning, is reproduced 
perfectly, with deviations of ~ 10 . On the other hand, the v e spectrum is not fitted very well. Without 
binning one can clearly see the steps caused by the fit. They disappear with binning, but the overall fit, 
particularly at higher energies, is not very good. That is caused by the choice of the fixed nuclear charge Z . 
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FIG. 5: (Color online) Average nuclear charge (Z) as a function of the (3 decay endpoint energy (in 0.25 MeV 
steps) for the 235 U fission by thermal neutrons. (For the definition of (Z) and its variance see text). The quadratic 
polynomial Z(Eo) that gave the best fit to the z/ e spectrum is shown by the dashed line. The lower panel shows, in 
log scale, the fission yield as a function of E en d po int- 



It turns out that different Z values contribute to different parts of the spectrum. This is shown in Fig. [5] 
There, for each bin with the endpoint Eq the quantity (Z) is calculated as 

(Z)(Eo) - ^fe4§¥§^ . A < Z >= «Z 2 > - m 112 ■ (3) 

The error bars in Fig. [5] are the dispersion A < Z > of the Z values within each bin. One can see clearly 
the slope of the (Z)(Eq) function; the (Z) values are decreasing with increasing Eq. At the same time, in 
the lower panel of Fig. [SJ one can see that the total fission yield decreases very rapidly with increasing Eq. 

The smooth dashed line in the upper panel of Fig. [5] shows the average nuclear charge Z giving fit to 
the v e spectrum that agrees with the original spectrum within 1% in the relevant interval of energies. It 
is represented by a quadratic polynomial in Eq as in some of the Refs. Q. The dashed line is above the 
corresponding < Z > due to the fact that in the eq.([T]) Z + 1 should be used. The quadratic polynomial 
fits of the form Z = a + b x Eq + c x Eq were used for all four nuclear fuels, resulting in deviation that 
does not exceed ~1% as shown in Fig. [BJ The corresponding parameters a,b,c are listed in Table I. The 
procedures (width of the slices and interpolating bins) used in the fits shown in Fig. [5] appear to be close to 
the optimum required for the precision conversion procedure. 

The calculated electron and D e spectra used for testing the accuracy of the conversion procedure were 
mutually consistent, but did not contain the small but nonnegligible QED and weak magnetism corrections. 
In conversion using the experimentally measured electron spectra these corrections should be included. The 
corresponding formulae can be found e.g. in Refs. 0,H, Q. 
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TABLE I: Parameters of the polynomial fit (Z) — a + b x Eo + c x Eq of the average nuclear charge as a function of 
the endpoint energy for the indicated isotopes 



nucleus 


a b c 


235-q 
239p u 

241 Pu 

238 tj 


50.0 0.825 -0.269 
50.0 0.825 -0.250 
50.0 0.825 -0.220 
50.5 0.825 -0.220 



0.02 



0.01 



-0.01 



-0.02 




3 4 5 6 

Antineutrino energy (MeV) 



FIG. 6: (Color online) Deviations of the v £ fits using the quadratic polynomial expression for nuclear charge depen- 
dence on the endpoint, Z(Eo) (see Table I), for the indicated nuclear fuels. Input spectra with 10 keV intervals were 
fitted in 100 keV wide slices which were then interpolated into 250 keV wide bins. 



IV. CONCLUSION 



Using the mutually consistent sets of the realistic electron and P e spectra corresponding to the fission of 
the nuclear fuels 235 U, 239 Pu, 241 Pu, and 238 U I have shown that the electron spectrum can be converted 
into the v e spectrum with an error that does not exceed ~1% in the energy interval 1-8 MeV. However, 
such accurate conversion can be obtained only if several conditions are met: 

• The slices into which the electron spectrum is divided are sufficiently fine. 

• The converted v e spectrum is smoothed out by binning in bins that are several times larger than the 
width of the original slices. 

• The optimum average nuclear charge (Z) is independently known as a function of the endpoint energy 
E . 

Not all these criteria were fulfilled in Refs. [f|, hence the larger uncertainties associated with the conversion 
procedure assigned there cannot be safely reduced. The electron spectra in [j| were published in 250 keV 
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wide bins. The number of bins and the number of slices for the conversion procedure were the same (or 
nearly so). The smooth electron spectra were expressed in the Fermi- Kurie representation which made the 
corresponding fits possible. And the v e spectra were not rebinned, thus the steps at the edges of slices 
were not smoothed out. However, as already pointed out above, the statistical and systematic errors of the 
measured electron spectra mostly exceeded the errors of the conversion procedure, so even if it would be 
somehow possible to minimize the uncertainty of the conversion procedure, it would affect the spectra of 
Refs. |5j only minimally. They remain the best existing source of the fission associated v e spectra. 

The main purpose of the present work is an examination of the conversion procedure per se. While there 
are no immediate plans to repeat and/or improve the results of the experiments in Refs. [a], there exist 



measurements of the (3 spectra of many individual neutron rich fission fragments (see, e.g. [10. These, 
in turn, can be converted into the v e spectra using the procedure described here, and combined with the 
fission yields to obtain the full v e spectrum. Similar measurements of individual short lived (and high Q 
value) fission fragments can be performed with the modern isotope separation facilities. Thus, significant 
improvements of the determination of the reactor D e spectra is possible. The present work is part of that 
effort. 
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